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Abstract. Discrete breathers are time-periodic, spatially localized solutions of the 
equations of motion for a system of classical degrees of freedom interacting on a lattice. 
We study the existence of energy thresholds for discrete breathers, i.e., the question 
whether, in a certain system, discrete breathers of arbitrarily low energy exist, or a 
threshold has to be overcome in order to excite a discrete breather. Breather energies 
are found to have a positive lower bound if the lattice dimension d is greater than 
or equal to a certain critical value d c , whereas no energy threshold is observed for 
d < d c . The critical dimension d c is system dependent and can be computed explicitly, 
taking on values between zero and infinity. Three classes of Hamiltonian systems are 
distinguished, being characterized by different mechanisms effecting the existence (or 
non-existence) of an energy threshold. 



PACS numbers: 45.05. +x, 63.20.Pw, 05.45.-a 
1. Introduction 

Discrete breathers are time-periodic, spatially localized solutions of the equations of 
motion for a system of classical degrees of freedom interacting on a lattice. They are 
also called intrinsically localized, in distinction to Anderson localization triggered by 
disorder. A necessary condition for their existence is the nonlinearity of the equations of 
motion of the system, and the existence of discrete breathers has been proved rigorously 
for some classes of systems [UEHHIIIIEIEIII!- In contrast to their analogues in continuous 
systems, the existence of discrete breathers is a generic phenomenon, which accounts for 
considerable interest in these objects from a physical point of view in the last decade. 
In fact, recent experiments could demonstrate the existence of discrete breathers in 
various real systems such as low- dimensional crystals jB], antiferromagnetic materials 9\, 
Josephson junction arrays [HUH], molecular crystals [12], coupled optical waveguides 
|13j . and micromechanical cantilever arrays [T3j . 

Properties of discrete breathers, as well as of some generalizations of discrete 
breathers, have been studied in detail in a large variety of different models. However, 
apart from existence proofs and studies of the spatial localization of discrete breathers 
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(which is typically exponential), only a few general results exist. Among these it is 
worth mentioning the remarkable result by Flach, Kladko, and MacKay ^5] on energy 
thresholds for discrete breathers in one-, two-, and three-dimensional lattices (and 
subsequent generalizations to systems with long range interactions [IS] and to partially 
isochronous potentials [T7j). Results on energy thresholds have practical relevance, as 
they can assist in choosing a proper energy range for the detection of discrete breathers 
in real experiments or computer experiments. It is found that, for Hamiltonian systems 
of infinite size, energy thresholds depend on the spatial dimension d of the system. A 
critical spatial dimension d c exists, such that 

• for a system whose spatial dimension d is smaller than d c , discrete breathers of 
arbitrarily low energy can be found, i.e., no energy threshold exist, 

• for a system where d > d c , an energy threshold exists, i.e., there is a positive lower 
bound on the energy of discrete breathers. 

The existence or absence of an energy threshold for discrete breathers beautifully 
explains certain observations made in discrete systems. One example are statistical 
properties characterizing the spontaneous formation of discrete breathers in cooled 
lattices [IB], a second example are the power spectra observed in thermalized lattices 
[T9] . In both cases, the presence of an energy threshold induces a qualitative change in 
the quantities under investigation. 

The result on energy thresholds in [T^] was derived under assumptions which 
turned out to be not general enough to cover all classes of systems of physical interest, 
and neither all of the different mechanisms leading to the (non) existence of an energy 
threshold. In this article, a much richer scenario is described, taking into account a larger 
class of systems and, most notably, pointing out three different mechanisms effecting 
the existence (or non-existence) of an energy threshold. The conditions leading to the 
occurrence of any of these mechanisms will be worked out, providing means to predict 
the existence of energy thresholds. A flavour of at least a part of these results has 
been given previously in a letter [201 - ^ n ^ ne present article, a more general situation is 
considered and details and derivations are given. 

There are three additional classes of systems, not covered in [T^j, for which energy 
thresholds and critical dimensions are discussed in the present article: Firstly, systems 
with nonanalytic potential are included in the analysis. Albeit rare, physical systems 
with such interactions exist and are under experimental investigation, for example in 
the field of granular media [21]. Secondly, systems in which no low amplitude discrete 
breathers exist are considered. Even for the class of systems with analytic potentials, 
these systems constitute a generic subclass, which accounts for their possible relevance 
in physics. Thirdly, results indicating the absence of an energy threshold in systems with 
no linear spectrum are presented. These kind of systems have attracted much interest 
recently, and strongly nonlinear tunable phononic crystals showing such a behaviour are 
currently being constructed [22]. 

The article is organized as follows: in section [2] some notation is fixed, defining the 
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classes of systems under consideration. Then, it is appropriate to distinguish between 
systems with and without a linear spectrum. Section |3] is devoted to the first category, 
and the existence of an energy threshold is traced back to the occurrence of a certain 
modulational instability. Systems without a linear spectrum are treated in section 0] 
Finally, in section 03 the results for the various cases are summarized, providing a 
compilation of the different mechanisms which effect the existence of an energy threshold 
and specifying the conditions which give rise to each of these mechanisms. 

2. Equations of motion 

We consider a hypercubic lattice Cd = {l,...,L} d in d 6 N spatial dimensions. For 
simplicity, we consider L G 2Z an even number. Each of the N = L d sites is labelled by 
a d-dimensional vector n G Cd, and to each site a state (p n , x n ) is assigned, where both, 
the momenta p n G M.f and the positions x n G W , are vectors of a finite number / of 
components. The dynamical properties of the system are gouverned by the Hamiltonian 
function 



with periodic boundary conditions in all spatial directions. V is called on-site potential, 
W 7^ interaction potential, and site n is assumed to interact with its neighbourhood 
Af n of nearest neighbouring sites on the lattice. 

The potentials V and W are both assumed to attain minima for zero argument, 
and, without loss of generality, V(0) = and W(0) = 0. In contrast to previous 
related work j2Sl EEH1 EI] , analyticity (in the sense of the existence of a Taylor series) 
of the potentials around their minima is not required. By removing this restriction of 
analyticity, additional models of physical interest are included, as for example granular 
media are described by means of nonanalytic interaction potentials [2"Tj . 

Several of the above restrictions on the system are imposed only for notational 
simplicity of what follows. Extending the computations to systems of different lattice 
geometries or interactions not only with nearest neighbours should be possible without 
any substantial change and are discussed to some extend in section 13.81 The case of 
long-range interactions is not included in the following analysis. A treatment of such 
systems can be found in ^B], and these results will be confronted with those for short- 
range interactions in section 

For the sake of readability, analytic calculations will be performed for the one- 
dimensional case 



again with periodic boundary conditions p n +N = Pn and x n+ N = x n , where momenta 
as well as positions are assumed to be one-component (/ = 1). The effect of the 
dimensionality on the results will be discussed in section 




(1) 




(2) 
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The existence of an energy threshold for discrete breathers will depend on the 
asymptotic behaviour of the potentials V and W close to their minima, and therefore 
we choose to represent V and W as power series in their arguments, 



\x\ r », 



X 



(3) 



with real (but not necessarily integer) exponents r M > 1. For convenience we consider the 
powers to be of increasing order, r M < r M+ i for all //. Some of the coefficients and/or 
M may be zero, and in this way the expansions in Q can also account for potentials V 
and W having different powers. Note that, by treating this class of potentials, we will be 
able to make statements about energy thresholds for the even larger class of potentials 
which have (j3J) as expansions around their minima. The potentials V and W in Q are 
chosen to be symmetric, and this restriction is made to keep the calculations simpler and 
the presentation more readable, as the number of cases to be distinguished is reduced 
in this way. For systems with linear spectrum and in the case of analytic potentials, 
the analysis of systems with asymmetric potentials can be found in [23]. These results, 
as well as the possibility of a generalization of our results to more general systems with 
asymmetric potentials, are briefly discussed in section 
The equations of motion following from (J2J) are 

x n + V'(x n ) + W'(x n - x n _i) - W\x n+1 - x n ) = 0, (4) 

where a dot denotes a total derivative with respect to time, a prime the derivative of 
a function with respect to its argument. It is convenient to introduce a transform to 
normal coordinates 

N 



N ^ 

n=l 



2ixl 



,1, i 



q 



5 6 fc i 2 ^ ••• 2 ; 



where the inverse transform is given by 



X n ^ G ^ Qq- 



(5) 



(6) 



The normal coordinates Q q (t), like the original ones x n (t), are functions of time, but 
this dependence will often be suppressed in the following. Rewriting the equations of 
motion (0J in terms of normal coordinates (0) yields 



Q q + F q (Q) = 0, 



where 



F q (Q) 



1 N 

N 2^ e 



n=l 



V 



-iq n 



and Q denotes a vector with entries Q q 



(7) 



+W ^(1 - <F)e-*»Q 9 J - W \^{ 



l)e~ lqn Q q 



(8) 
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3. Systems with linear spectrum 

In this section, we want to refer to systems with interaction potential W(x) = 
4r£ 2 + 0(x s ) quadratic in leading order, i.e., (f>o ^ and s > 2. (O denotes Landau's 
order symbol.) The on-site potential, in contrast, can be either zero (V = 0), or of 
quadratic or higher order, V(x) = 0(x r ) with r > 2. In the notation of the expansions 
in (J3J) this means r = 2 and 0o 7^ 0. This setting of parameters allows linearization of 
the equations of motion (jlj) with a non-trivial result, 

Q q + u*Q q = 0, (9) 

where 

u\ = ^o + 40o sin 2 (|) (10) 

denotes the frequency of the linear mode q. The set of all the u q is called the linear 
spectrum. The differential equations © obviously decouple, having simple sine waves as 
solutions. In the following we will be interested in solutions of the nonlinear equations of 
motion (j3J) in the limit of small oscillation amplitudes in which the linearized equations 
of motion Q are approached. 



3.1. Main ideas of this section 

For generic Hamiltonian systems, discrete breathers, like all periodic orbits, occur in 
one-parameter families. Some typical choices of quantities to index such a family are 
the energy of a discrete breather, its frequency, or its amplitude measured at the site of 
maximum amplitude. In numerical studies of discrete breathers in a variety of systems 
with linear spectrum, the following scenario has been observed: consider some discrete 
breather with frequency u DB outside the linear spectrum. When following the family of 
discrete breathers and approaching the edge u; edge of the band formed by the frequencies 
of the linear spectrum, it may or may not be the case that the breather amplitude 
A goes to zero. The reverse conclusion, however, can be drawn: if discrete breathers 
of arbitrarily low amplitude exist, they cannot be found elsewhere but in the limit 
u DB — > u; odgc [T3]. Furthermore, in this case where 

lim A db (cj db )=0, (11) 

it is observed numerically that in this limit the localization of the discrete breather 
becomes weaker and weaker. This fact is illustrated in figure ^ where, for the example 
of the a Fermi-Pasta-Ulam chain and discrete breathers of various frequencies, the 
maximum values of the amplitude are plotted for a number of lattice sites. More 
precisely, the localization strength approaches zero in the limit u UB — > co> odgc , and a 
transformation of a discrete breather into a spatially extended solution (plane wave) 
appears to take place. This observation gives rise to the following 

Hypothesis 1 Consider a system with linear spectrum as specified above. Assume that 
in this system there exist low amplitude breathers, i.e., a family of discrete breathers for 
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Figure 1. Shape profiles (maximum values of the amplitude) of discrete breathers, 
centered at lattice site 0, in a Fermi-Pasta-Ulam chain of N = 99 degrees of freedom 
(the equations of motion are given in (|48|l with parameter value a = 0). For a system 
which — like this one — has a linear spectrum, the localization is exponential (straight 
lines in the logarithmic plot). When the maximum breather amplitude is lowered (from 
+ to ■), the localization strength decreases and the solution attains more and more 
the shape of a plane wave. Lines are merely drawn to guide the eye. 

which (G2P holds. Then low amplitude breathers emerge from band edge plane waves by 
means of a tangent bifurcation. 

A tangent bifurcation by definition takes place when a pair of Floquet multipliers 
of the linear stability analysis of the periodic orbit happens to collide at +1 on the 
unit circle. The periodic orbit emerging from the bifurcation appears to be of the same 
frequency as the one it originates from, and this is exactly the situation we observe for 
discrete breathers stemming from band edge modes. 

The following three results or observations may serve to substantiate this hypothe- 
sis. Firstly, in a recent existence proof of discrete breathers due to James P , it is shown 
for a large class of one-dimensional systems that discrete breathers of small amplitude 
exist and that they emerge from a bifurcation of a band edge plane wave. Secondly, 
Flach has shown in that the periodic orbits emerging from the tangent bifurcation 
lack a certain permutation symmetry. This is consistent with the periodic orbits being 
discrete breathers, but inconsistent with plane wave-type solutions. Thirdly, to the 
knowledge of the author, in the numerous numerical studies of discrete breathers, there 
has never been observed a family of discrete breathers in a system with linear spectrum 
which, in the limit (if it exists) of the amplitude going to zero, is not subject to the 
attenuation of the localization strength as described above. 

The above hypothesis allows to tackle the problem of the existence of an energy 
threshold for discrete breathers in the following way (which is the outline of the rest 
of this section): instead of directly treating discrete breathers, the analytically more 
easily accessible band edge plane waves are considered. Since it is the behaviour at low 
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amplitude which will be of importance, this can be done by employing perturbation 
techniques. Performing a Floquet analysis of these band edge plane waves, criteria for 
the appearance of a tangent bifurcation are found and the corresponding bifurcation 
energy is obtained. Then, the above hypothesis allows to relate this bifurcation energy 
to the energies of discrete breathers. In case the bifurcation energy is vanishing, no 
energy threshold is present. If the bifurcation energy is positive, an energy threshold 
has to be overcome in order to excite a discrete breather. 

Note that there exist systems in which no low amplitude breathers are present 
and consequently the above hypothesis does not apply. This special case is deferred to 
section E33 whereas for the following sections we assume (fTTj) to be valid. 



3.2. Band edge modes 

At the band edges of the linear spectrum, there exist two particularly simple periodic 
orbits Q 1 and Q u , each with only one mode excited, 

I: Q q = Vg^O and Qo+V'{Q ) =0, 
II: Q q = Vg^vr and Q w +V'(Q w )+2W'(2Q n ) = 0, 

where the excited mode is determined as the solution of the respective differential 
equation. Q 1 is called the in-phase mode, characterized by identical oscillations at 
all sites, i.e., x n+ i(t) = x n (t) for all n 6 {1, N} at all times t. As is obvious from the 
differential equation for Q in (fl^jll. the presence of an on-site potential V is essential 
for the existence of an in-phase mode. Q 11 we will refer to as the out-of-phase mode, 
where each two neighbouring sites oscillate in opposition of phase, i.e., x n+ i(t) = —x n (t) 
for all n G {1,...,N} at all times t. Exemplarily, the following computations will be 
presented for the in-phase mode only. 

Of a uniformly valid perturbative solution for the differential equation (jl2|) I. for 
our purposes the leading term in the energy e is sufficient, 

Qo(t) = J—coBp(X)t + <p} + h.o.t, (13) 
V 



which can be obtained for example by means of a Lindstedt-Poincare expansion [23j at 
low energy 

e(Q 1 ) = lQl + V(Qo). (14) 

<p is an arbitrary phase shift which will be set to zero in the following, and higher order 
terms (h.o.t.) in the energy have been neglected. To render the expansion uniform, the 
nonlinear oscillation frequency 



1+ *V i 9 X + h.o.t. 



2 



(15) 



has to be determined up to first order in 

2e \ (n-2)/2 



*=(-] (16) 
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where T is the Euler gamma function and higher order terms in X have been neglected. 
This series expansion can be derived either by standard perturbation techniques [21] or 
using a method put forward in [T7j . 

3.3. Tangent bifurcation 

Our aim is to determine the energy at which a band edge mode undergoes a tangent 
bifurcation. Such a bifurcation implies a pair of Floquet multipliers of the linear stability 
analysis of the periodic orbit to collide at +1 on the unit circle, and is hence accompanied 
by a change in stability. In order to investigate the stability properties of the band edge 
modes, a small perturbation is considered by substituting 



in the equations of motion 0. In the case of the periodic orbit Q 1 , this substitution 
gives rise to equations of motions for the perturbations S q which, up to linear order in 
5 g , read 



a periodic function in time. This type of differential equation (ordinary, homogeneous, 
second order, with periodic coefficient) is termed Hill equation, and a number of 
theorems concerning stability are at hand. Let us briefly recall some properties of 
the Hill equation (for details see |2Sj) arid outline how they can be exploited in order to 
determine the bifurcation energy. 

It follows from Floquet theory that, depending on the values of the parameters 
present in (fT£J) . the Hill equation has either stable, unstable, or periodic solutions. 
Considering a suitable two-parameter space — in our case we choose the linear frequency 
u q and the energy e, where the latter in turn determines V" — the regions of stable 
solutions display a particular shape in parameter space, the so-called Arnold tongues 
(see figure 11-3 in [21] f° r an illustration). Regions of stable solutions are separated 
from regions of unstable solutions by lines of periodic solutions. Hence, the parameter 
values at which the band edge mode Q 1 bifurcates can be obtained by determining the 
values at which any of the 5 q has a periodic solution. The tangent bifurcation we are 
looking for appears at the smallest non-zero energy for which a periodic solution (|T3*j) 
of orbit I and a periodic solution of the corresponding variational equation (|TH|) . both 
with the same period, exist. In the following, by determining an energy expansion of the 
frequency f2(e) of the periodic solution of the variational equations (fTHj). and equaling 
it with the frequency Q(e) of the periodic orbit ()15j) . the bifurcation energy will be 
obtained to leading order. 



Qq — > Qq + 0~q 



(17) 



5 q (t)+[V"(t)+u; 2 g }5 q (t) = 



(18) 



where 
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The periodic solutions, as well as the corresponding parameter values, can be 
obtained by perturbation methods. In order to apply the Lindstedt-Poincare technique 
[21], we introduce the transformation r = Q(X)t into QlKj) . yielding 

n 2 (X)5 q (X, r) + [V"(X, r) + uj 2 q ] 5 q (X, r) = 0, (20) 

where a dot now denotes a total derivative with respect to r. This change of variable 
transforms V" in a function which is ^-periodic in r. The energy dependence has been 
noted explicitly by the variable X. Inserting the Ansatze 

5 q (X, r) = 6 qfi (r) + 5 qA (r)X + h.o.t., (21) 

n 2 (x) = ng + n?x + h.o.t., (22) 

as well as the series expansions (fTH|) and (|TK|) in the Hill equation (]2"Uj) . terms of equal 
order in X can be collected to obtain the differential equations 

order X° : fig^o(r) + ^\, (t) = 0, (23) 
order X 1 : fljfo,i(r) + u 2 q 5 q}1 (r) + O?5 9 , (r) + Vf (r)(J,, (r) = 0, (24) 

where the 7r-periodic function 

V'((t) =w 1 (r 1 -l)|cosr| n - 2 (25) 

has been defined. The solution of the zeroth order equation (}2"3*|) is 

S qJ0 (r) = A cos + B sin , (26) 

where A and 5 are complex valued constants. It follows from Floquet theory that, in 
order to obtain periodic solutions of period 2ir of the first order equation (|2*4"|) in the low 
energy limit, it is sufficient to consider only the first two terms of a Fourier series of the 
driving term, 

oo 

Vf(r) = J2 C kCos{2kr) (27) 

k=0 

with coefficients 

2 Vl r(**±i) 

°° = 7^T#T' (28) 



2 i;i(ri-2)r(n) 

Ci - 2 r(^)r(t) • (29) 

Inserting the truncated Fourier series as well as the zeroth order solution ([26)1 . the first 
order equation ([2*4*jl reads 



fig(y,,i(r) + ^ 2 <W(t) + [£> + Ci cos(2r)] 5 q , (r) = 0, (30) 

where 

D = q,-|Spy. (3D 
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The differential equation (J3U|) is linear, and its solutions, which are in general aperiodic, 
can be obtained. Then, the periodic ones are found by eliminating the secular terms in 
the solution. This elimination is achieved only under certain conditions on the frequency, 
where the ones corresponding to the tangent bifurcation we are looking for are given by 

Ql = u 2 , (32) 

= Co + y- (33) 

Equalling the truncated series expansions of the squared frequencies Q 2 and Q 2 , and 
solving for the variable X, the critical values 

/ 2\ 2 /(n- 2 ) 

<«=i{^cr) < 34 > 

of the energy at which a tangent bifurcation of orbit I may occur are obtained. Note that, 
in the preceding derivation of the bifurcation energy, we have simultaneously considered 
N uncoupled equations (fTH|) for the perturbations S q , labelled by the index q. From <5 
no bifurcation can arise, as this component simply effects a time shift of the periodic 
orbit |2Sj- A change of stability in any of the remaining equations, however, may give 
rise to a bifurcation. We are interested in the tangent bifurcation giving the lowest value 
of the bifurcation energy (|H4jl . which implies the choice q = 2ir/N . For a large number 
of lattice sites N, we can write the corresponding linear frequency (fTU|) as 

</n = ^ + 47r 2 o iV- 2 + 0(N~ 4 ). (35) 

Inserting this expansion as well as equation (|29j). the bifurcation energy of orbit I reads 



<2,r/iV - 2v 



7T 2 o r 1 r(n-2 ) r(f)^ /(ri - 2) 



(36) 



-N 2 ViT(n) 

in leading order of iV -1 . Since Vq, <fto > and since we have to demand a positive 
bifurcation energy, the tangent bifurcation can only occur if 

V! < 0. (37) 

A similar analysis for the out-of-phase mode (periodic orbit II) yields the expression 

2/(n-2) 

(3* 



n Vq + 40o 



^5/2 0orir(I1 . 



2N 2 ( Vl + 2'i0 1 )r( , 

for the bifurcation energy. We have v o > and <f) > 0, leading to the necessary condition 
v x + 2 ri 0i > (39) 

for a tangent bifurcation to take place. Similar results in |2Sj, derived for the case of 
analytic potentials V and W, are recovered by setting r\ = 4. It was already noted 
in the same reference that conditions ()37)1 and (|3T)|) mirror the fact that, for a tangent 
bifurcation to take place, it is necessary that the frequency of the band edge mode is 
repelled from the linear spectrum of the system with increasing energy. 
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3.4- Higher spatial dimensions 

The computations and results of the preceding section can be generalized to systems of 
higher spatial dimensions d. The changes in the expressions of the bifurcation energies 
(JUS)) and (|3~8J) are of small extent, and for our purposes it will only be important to 
note how the proportionality between the bifurcation energies e c and the system size N 
changes. 

For arbitrary spatial dimension d, the normal coordinates (JHJ) are labelled by a d 
dimensional vector q = (q 1 ,...,q d ) with = an d k E {—L/2 + 1, ...,L/2} for all 
i. Then, choosing the value of q giving rise to the lowest bifurcation energy, a factor 
£4/(2-n) j g obtained in the expressions of the bifurcation energy instead of jV 4 /( 2 ~ r i). 
Since iV = L d , this is incorporated into (pffij) and by replacing N by N l l d . In this 
way for both, the in-phase as well as the out-of-phase mode, the proportionality 

e c txiV 4 /W 2 -^)], (40) 
connecting the bifurcation energy e c with the system size N, is obtained. 



3.5. Energy thresholds for discrete breathers 

The hypothesis of section I3~T1 allows us to relate the bifurcation energy of the band edge 
plane wave to the minimum energy accessible for a discrete breather. It follows from 
and ((Hj) that the bifurcation energy determined above is an intensive quantity, i.e., on 
the scale of total energy per particle, and we will now turn to its extensive counterpart 

E c = Ne c . (41) 

Considering the proportionality (j4*Uj) in the limit of infinite system size N — > oo, it is 
found that 

lim E c oc lim N 1 ^^ { = ° * ^ (42) 
JV— oo N^oo I >0 if d>d c , 

with critical dimension 

dc = (43) 

Hence, for systems of spatial dimension d smaller than the critical value d c , discrete 
breathers of arbitrarily small energy can be found, whereas for systems where d > d c 
an energy threshold has to be overcome in order to excite a discrete breather. Equation 
(I43J1 is in perfect agreement with the critical dimension d c — 2 obtained in [15J for 
(generic) analytic potentials with a non-zero quadratic and quartic term, which, with 
our restriction to symmetric potentials, corresponds to a nonlinearity of leading order 
T-i = 4. 

This result for the critical dimension can be recovered by considering a continuum 
approximation of the discrete system as explained in pH]- In this approximation, the 
breather energy Ef, can expressed as an integral (equation (6) of reference [T5] ) which, 
in the limit of small breather amplitudes A, yields the proportionality Eg oc A^~ zd ^ 2 . 
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For our choice of symmetric potentials V and W in (jHJ), the detuning exponent z (see 
[To] for a definition) is equal to T\ — 2. Then, considering in the limit A — ■> 0, the 
expression (|43j) for the critical dimension cf c , distinguishing between finite and diverging 
limiting values for E b , is reproduced. 

3.6. Intuitive understanding 

The mechanism which gives rise to an energy threshold for discrete breathers in high 
enough spatial dimension can be understood on an intuitive level. When lowering the 
amplitude of a discrete breather towards zero, two different mechanisms, competing with 
respect to their effect on the breather energy, take place: the localization strength of the 
discrete breather tends to become weaker, leading to an increase in energy, while the 
mere lowering in amplitude causes a decrease in energy (see figure ^ for an illustration). 
Depending on the respective strengths of these effects, an energy threshold may or may 
not exist. In the limit of weak localization, the degrees of freedom far from the breather's 
centre gain in importance. Their contribution to the energy depends on their 'number', 
which in turn depends on the spatial dimension of the system, and it is in this way that 
d enters the game. 

3. 7. Numerical confirmation of the results 

Numerical computations have been performed in order to confirm the above result on 
energy thresholds for discrete breathers in systems with linear spectrum. By numerical 
continuation of periodic orbits from an anti-continuum limit discrete breathers on 
finite lattices can be computed numerically up to machine precision. In particular, it is 
possible to continue a discrete breather along its family while varying a parameter. In 
doing so for a set of frequencies approaching the band edge of the linear spectrum, the 
dependence of the breather energy on its amplitude, measured at the site of maximum 
amplitude, can be determined. In order to test our main result, the existence or non- 
existence of an energy threshold of discrete breathers under certain conditions, we will 
confront such data for two exemplary systems, one displaying an energy threshold, the 
other not. 

Consider a two-dimensional system (d = 2) of Fermi-Pasta-Ulam type, i.e., with 
zero on-site potential V = 0. The interaction potential W in the Hamiltonian ((TJ) is 
chosen as 

W(x) = -x 2 + —\x\ ri , (44) 
2 ri 

giving rise to equations of motion of the form 

X n -\- ^ ^ (•£»! ^m) [l \%n •Em\ ] 0, (45) 

where Af n denotes the set of nearest neighbours of lattice site n. Then, following ()42j) 
and (|43|). no energy threshold is present for r\ = 3, whereas for the choice of r\ = 4 
a threshold is obtained. These predictions are confirmed by the numerical results for 
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Figure 2. Breather energy versus amplitude for two-dimensional Fermi-Pasta-Ulam 
systems <|45|) with periodic boundary conditions for the cases r± = 3 and r\ = 4. A 
lower bound on the breather energy is observed. For r\ — 3, with increasing system 
size N this bound converges towards zero, whereas for r\ = 4 it approaches a non-zero 
value. 

finite systems with periodic boundary conditions presented in figure El For any finite 
number of lattice sites, a lower bound on the breather energy is observed. For r\ = 3, 
with increasing system size N this bound converges towards zero, whereas for r\ = 4 
it approaches a non-zero value. Note that, for the system sizes considered, the range 
of validity of the large N approximation (}3*H|) has not yet been reached, as, in this 
approximation, no system size dependence of the energy threshold is expected in the 
case of ri = 4. 

3.8. Generalizing to larger classes of systems? 

The above result on energy thresholds for discrete breathers has been obtained for the 
class of systems defined in section |21 However, several of the restrictions have been 
imposed merely to keep the presentation clearer and for notational convenience, and we 
will discuss in the following which ones presumably could be removed. 
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It is obvious from (|42j) that the existence of an energy threshold is ruled by the power 
of the system size N in the expressions of the bifurcation energy (j3T?j) and ([38)1 . This 
power, in turn, can be traced back to the powers of the energy X in the series expansions 
of the frequency (fT3j) of the band edge mode, and of the frequency (jUJ) at which 
parametric resonance occurs. Although we have not explicitly done the calculations, 
we do not expect the powers in these expansions — and hence the energy thresholds for 
discrete breathers — to vary under a number of generalizations like 

• other than hypercubic lattice geometries, 

• states (p n , x n ) with more than two components (/ > 1 in the notation of section EJ), 

• interaction with a finite neighbourhood larger than nearest neighbours only, 

Energy thresholds for discrete breathers in systems with Hamiltonian functions of other 
than the standard form (JI}, like the discrete nonlinear Schrodinger equation, have been 
considered in [TH] and, on a rigorous level (i.e., without invoking our hypothesis 1), 
in [21]. Systems with asymmetric V and W are discussed in for the case of 
analytic potentials. There it is found that allowing for asymmetric contributions in 
the series expansions (JHJ) changes the expression for the bifurcation energy E c , but it 
does not change the critical dimension d c , and we would assume the same to happen for 
nonanalytic potentials. What can be modified by the inclusion of an asymmetric term, 
however, is whether low amplitude breathers do or do not exist in a certain system, and 
this issue is discussed in the next section. 

For systems with long range interactions, the existence of energy thresholds is 
investigated in 16J. There it is found that the presence of long-range interactions 
enhances the appearance of an energy threshold, as shown for one-dimensional systems 
with analytic potentials, in which, in contrast to the short-range case, an energy 
threshold is observed. 

3.9. Systems without low amplitude breathers 

The above analysis leading to our result on energy thresholds for discrete breathers 
is heavily relying on the hypothesis stated in section where it is assumed that 
breathers of arbitrarily low amplitude A DB exist. This condition is fulfilled in many 
cases, but exceptions not only do exist, but they even appear to be generic. The class 
of systems without low amplitude breathers is defined by the property 



for all families of discrete breather present in the system when the breather frequency 
u DB approaches the frequency ou cdga of any of the band edge modes. First we want to give 
an example of a system without low amplitude breathers and study its characteristics. 
As a second step, some speculations are made on the general conditions leading to the 
absence of low amplitude breathers. 



lim 




(46) 
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Consider a Fermi-Pasta-Ulam (FPU) chain, defined by Hamiltonian (J2J) with 

V = and W(x) = ^x 2 + jjx 3 + ^x\ (47) 
which gives rise to the equations of motion 

i n + [(x n - x m ) + a(x n - x m ) 2 + (x n - x m f] = 0. (48) 

m=n±l 

For this system, two rigorous proofs of the existence of discrete breathers have been 
published. The first one is due to Aubry, Kopidakis and Kadelburg using a 
variational method. These authors proof existence of discrete breathers for any strictly 
convex interaction potential W, which corresponds to a parameter a with \a\ < v3 in 
(jUJ). A second proof, applying a method of centre manifold reduction, was obtained 
by James [Sj, where existence of discrete breathers has been shown for values of a 
obeying \a\ < |v3. Furthermore, as this latter method guarantees to catch all small 
amplitude solutions present in the system, the existence of small amplitude breathers can 
be excluded for \a\ > |a/3. Recapitulating, for the FPU chain (}4*%j) and parameter values 
|\/3 < a < Vs, existence of discrete breathers can be proofed, but small amplitude 
breathers in the sense of (jllj) are clearly absent. Knowing that low amplitude breathers 
do not exist, it is a triviality to infer the existence of an energy threshold for discrete 
breathers. 

To illustrate this result, numerical data, obtained by the continuation method 
described in section l3~7| are presented. The two cases a = 0.4 < ~v3 and |\/3 < a = 
1.2 < v^3 are considered, which allows to confront the behaviour of systems with and 
without low amplitude breathers. Varying the breather frequency uj ob and approaching 
the band edge of the linear spectrum, the breather amplitude at the site of maximum 
amplitude is monitored. The result is presented in figure El showing the amplitude to 
converge towards zero for u DB — > 2 in the first case (a = 0.4), whereas a non-zero value 
is approached in this limit for the second case (a = 1.2). 

Some speculations can be made regarding the general conditions under which 
exclusively large amplitude breathers (J46j) are observed, implying the existence of an 
energy threshold for discrete breathers. At the end of section 13.31 we noticed already 
that, for a tangent bifurcation of the band edge mode to take place in the class of systems 
considered, the expansion coefficients of the potentials Q have to fulfill the inequalities 
(j3Tj) and (J3TJJ) for the in-phase and the out-of-phase mode, respectively. For the FPU 
system (HU), the corresponding inequality was already obtained in |2SI, equation (3.24), 
yielding a < |a/3 in our notation. This inequality obviously coincides with the one 
deduced from the rigorous results flj and jH], identifying the values of a providing low 
amplitude breathers. From this observation, one might deduce the following 

Hypothesis 2 //, in a system with linear spectrum, no tangent bifurcation of any of 
the band edge modes take place, then discrete breathers of low amplitude in the sense of 
im\) do not exist. 
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Figure 3. Breather amplitude versus frequency for system l|48[) with N — 99 lattice 
sites, free boundary conditions, and frequencies close to the band edge of the linear 
spectrum (which is at frequency 2.0). With parameter value \a\ = 0.4 < ^V3, 
breathers of arbitrarily low amplitude are observed, whereas for < \a\ — 1.2 < ^/i 
the breather amplitude, and hence the breather energy, does not tend to zero when 
approaching the band edge of the linear spectrum. 



This hypothesis, if correct, allows to deduce the existence of an energy threshold 
for discrete breathers from the mere absence of a tangent bifurcation of the band edge 
modes. 

4. Systems with no linear spectrum 

In analogy to the definition at the beginning of section El by the expression "systems 
with no linear spectrum" we want to refer to systems whose potential functions V and/or 
W of the Hamiltonian are not quadratic in leading order. In the notation of the 
expansions in (JHJ) this means simply r ^ 2. The peculiar properties of such systems 
have attracted much interest recently, and strongly nonlinear tunable phononic crystals 
showing such a behaviour are currently being constructed [22 3 . 

In contrast to systems with linear spectrum where discrete breathers are known to 
be exponentially localized in space pUEE], a localization even stronger than exponential 
is found in systems with no linear spectrum [SEIESIED]- An upper bound on the breather 
amplitude is given by the inequality 

\xi\<aexp(-\i\/b)c {dtm (49) 

for some a,b,c,d G R in one-dimensional systems jHOj, where the breather is centred 
at site i — 0. However, arguments by which these results for the spatial decay can 
be obtained should be equally applicable in higher spatial dimension, and numerical 
computations confirm this reasoning, finding superexponential localization in higher 
dimensional systems. 
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Again, as in the preceding sections, we are interested in the low amplitude behaviour 
(if existing) of discrete breathers, which in turn affects the existence of an energy 
threshold. The mechanism as described in section 01 the emergence of discrete breathers 
from a bifurcation of a band edge plane wave, cannot take place here simply due to the 
absence of a linear spectrum. Nevertheless, it might in principle be the case that, when 
lowering the amplitude of a discrete breather in a system with no linear spectrum, 
the localization strength decreases, giving rise to an energy threshold in some spatial 
dimension. For lack of any better idea, some numerics was performed, suggesting for 
the (few) cases considered that 

(i) low amplitude breathers appear to exist, and 

(ii) the localization strength of discrete breathers remains constant in the low amplitude 
limit. 

As an example, data are presented for a Klein-Gordon type system with Hamiltonian 
(JTJ and potentials 

V(x) = -x 2 and W{x) — -\x 
giving rise to the equations of motion 

X n X n -\~ ^ {x m *^n)|*^m *£n|; 

where Af n denotes the set of nearest neighbours of lattice site n. Plotting the shape 
profiles of various discrete breathers obtained numerically for a one-dimensional system 
in figure EJ the above points (0) and (0) are corroborated. Inspired by the numerical 
observations, we formulate the following 

Hypothesis 3 Consider a system with no linear spectrum as specified above which 
supports discrete breathers. Then, low amplitude discrete breathers in the sense of (Q7J) 
can be found. Furthermore, the localization strength is approximately constant within a 
family of discrete breathers. 

Note that this hypothesis is based only on numerical observations in a few examples 
of systems, both of FPU and Klein-Gordon type. From this hypothesis, the absence of 
an energy threshold for discrete breathers in arbitrary spatial dimension is deduced 
immediately. For illustration, numerical data are presented in figure El for system (|51j) 
in spatial dimensions d = 1,2,3,4. Absence of an energy threshold is observed in any 
of these dimensions already for small finite systems. This is clearly distinct from the 
behaviour observed in figure 121 for systems with a linear spectrum where, if at all, the 
threshold vanishes in the limit of large system size only. 

5. Summary of the results 

We have studied the existence of energy thresholds for discrete breathers in a large class 
of Hamiltonian systems. Breather energies are found to have a positive lower bound if 



(50) 
(51) 
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Figure 4. Shape profiles (maximum values of the amplitude) of discrete breathers, 
centered at site 0, in a system with no linear spectrum, namely the one-dimensional 
version of a Klein-Gordon type system with equations of motion l|51|) and N — 99 
degrees of freedom. The localization strength is superexponential and remains 
unchanged when varying the maximum breather amplitude (or, equivalently, the 
breather frequency w DB - Note the twofold logarithmic scale of the ordinate 
(— ln[— ln(-)]). Lines are merely drawn to guide the eye. 
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Figure 5. Energy versus amplitude for discrete breathers in a system with no 
linear spectrum. Independently of the spatial dimension d — 1,2,3,4 (from top to 
bottom), arbitrarily low breather energies are observed. The data were obtained for 
systems consisting of N = 8 d lattice sites, but, due to the strong (superexponential) 
localization, are indistinguishable on this scale from those of larger systems. 



the lattice dimension d is greater than or equal to a certain critical value d Cl whereas no 
energy threshold is observed for d < d c . The value of d c depends on the class of systems 
under consideration, and we can distinguish the following three cases: 

(i) Systems with continuous linear spectrum which support low amplitude breathers: 
discrete breathers of low amplitude emerge from a tangent bifurcation of a band 
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edge mode. The bifurcation energy can be calculated explicitly and is related to 
a bound on the breather energy. For all finite systems, an energy threshold for 
discrete breathers is observed. For systems of large sizes N, this energy threshold 
scales as jV 1-4 /^ 1-2 )], resulting in a critical dimension d c = 4/(r x — 2) for infinite 
systems. (See section EJfor a definition of ri.) 

(ii) Systems with linear spectrum which do not support low amplitude breathers: no 
tangent bifurcation of a band edge mode is observed, and apparently no other 
mechanism for the creation of low amplitude breathers exist. As an immediate 
consequence of the absence of low amplitude breathers, an energy threshold for 
discrete breathers exists for finite as well as infinite systems. 

(iii) Systems with no linear spectrum: the localization strength of discrete breathers 
does not vary significantly with the frequency, and strongly localized low amplitude 
breathers appear to exist. This brings about the absence of an energy threshold for 
discrete breathers for finite as well as infinite systems, i.e., d c = oo. 

Note that the above results are derived on the basis of hypotheses 1 to 3 formulated in 
sections |3] and 0] The hypotheses are corroborated in part by analytical results, in part 
by numerical observations. 
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